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Abstract 

In previous work, we have argued that measurement using a radar can be viewed as taking the 
expected value of an operator. The operator usually represents some aspect of the characteristics of 
the object being tracked (such as Doppler, distance, shape, polarization, etc.) that is measured by 
the radar while the expectation is taken with respect to an optimal matched filter design process 
based on the waveform broadcast by the radar and a receiver which is optimized to a specific 
characteristic of the object being tracked. With digital technology, it is possible to produce designer 
waveforms both to broadcast and to mix with the return signal, so it is possible to determine the 
maximum of the expectation of the operator by proper choice of the received signal. We illustrate 
a method for selecting the choice of the return signal to detect different "target operators" using 
perturbation theory based on the Matched Filter Principle and illustrate it with different operators 
and waveforms. 
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In t 

Radar" 20|, Woodward introduced the ambiguity function as the means to solve the mea- 
surement problem of radar. The measurement problem of an active sensor is to design a 
waveform to be broadcast by a radar or sonar, to maximize the receiver response to the 
signal which has interacted with an object. The solution proposed by North[13| during 
World War II is the "matched filter", which correlates a known signal template with what 
is received in a return signal to detect the presence or absence of the template in the un- 
known received signal. This is exactly equivalent to convolving the unknown signal with 
the complex conjugate of the time-reversed version of the known signal template; this is 
called cross-correlation. Therefore, as has been shown in many texts[19|, the matched filter 
is the optimal linear filter for maximizing the signal to noise ratio (SNR) in the presence of 
additive noise. 

In radar or sonar, a known signal is sent out and the reflected signal from the object 
(which is a function of the distance to the object, the relative speed of the object and the 
broadcast frequency of the radar), can be examined at the radar receiver for the common 
elements of the out-going signal in the return signal, which, when optimized is a multi- 
dimensional matched filter or ambiguity function. The broadband form the return signal is 
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s(at — r), where, r = — , is the delay 

C-V R 1-p 

Here c is the speed of propagation and vr is the radial velocity of the object. 

There are two forms for the ambiguity function, the more general form is the wideband 
(WB): where the return signal can be modeled as a delay in time of the broadcast signal. 
The wideband ambiguity function, Xwbj which has the return signal modeled as both a 
dilation and delay of the broadcast signal 

/oo 
e- it0J s*(t)s N , WB (t)dt, (2) 
■oo 

where s*(t) means complex conjugate of the broadcast signal. The ambiguity function is 
used to design radar signals so that they have desirable properties useful for various kinds of 
radars (Leavon is a current up to date resource 12]). We propose a way to think about the 
ambiguity function which is different than the way Woodward presented it. This approach 
suggests the ambiguity function can be thought of as the expectation value of an operator 
that is connected to the delay and dilation properties associated with the Doppler effect [10 ]. 
Thus, the sensor measurement problem can be cast in a more abstract setting, which treats 
interaction between the waveform and the target as an operator acting on the waveform. This 
approach can be termed the operator approach and it can be viewed as an abstraction of the 
quantum mechanical formalism applied to a classical setting. This approach underlies the 
time- frequency approach to signal processing that has been championed by Cohen jj]. Using 
this approach, we examine the operator viewpoint for both single and multi-dimensional 
operators acting on a signal by the interaction process. In particular, we propose that the 
cross-ambiguity function for certain operators can be used to amplify the return signals. 
We illustrate this for several operators, show under what conditions this amplification can 
occur, and discuss how the cross-amplification signal can be constructed given knowledge of 
the interaction operator and the broadcast signal. The result of this approach is to suggest 
a way for recasting problems in signal processing when we have sufficient knowledge of the 
interaction of the broadcast signal. 



3 



II. OPERATOR APPROACH 



The notation for the inner product of two signals, r(t) and s(t) that is used throughout 
the paper is 

POO 

(r(t),s(t))= / r*(t)s(t)dt, (3) 



while the Fourier transform, J 7 , of a signal s(t) is[16| 

/oo 
e~ itu s(t) dt = (e ituJ ,s(t)) , (4) 
-oo 

and the inverse Fourier transform, is 

s(t) = F- l S(u) = — / e ltuJ S(u)duj = < e~ ituJ , S(u)) . 

A function of time which is translated by amount, r, can be written as (using the Taylor 
expansion of function D = 4| 



(if ; 



s (t + t) = e rD s{t) = e i y- lT -t) s{t) = e lW s{t) (5) 
The form of the narrow band (N) ambiguity function Xn( u i r )> can be recast as 

Xn (lu, t)= f e~ ltw s*(t)s(t -r)dt= (s(t)e itul , e~ iT ™ s(t)\ . (6) 

J — oo 

From the Doppler effect perspective, translation is the operation of the frequency operator 
on the signal, where r is the total distance a signal travels to an object, is reflected, and 
then returns to the receiver. The expected value associated with observable, A, for a signal 
s (t) is 

(A^ = jA\s{t)\ 2 dt = J s*(t) As(t) dt=(s(t),As(t)y (7) 
Thus, the narrow band ambiguity function can be written using this definition as 



«(*) 

We can thus interpret e ±trW as a translation operation acting on function s(t) which moves 
the time t — > t±r. This way of considering measurement in radar is a natural continuation of 



the viewpoint that started with Gabour[8| and extended by Woodward [20| and Vaidman[18l 
for considering measurement in radar. 
The time operator, T, is 

i du 



while the frequency operator is 

>V=~. (9) 

i dt w 

It is understood that these operators act on signals and that 

#"'W = (i£l" •<■*>■ (10) 

A very useful calculation trick is based on a modification of Parceval's theorem for an 
unnormalized signal: 

£=<*(*),*(*)> = <1>. W 

= ^ F ((5(a/)e iw 'M),<5( W )e <wt > l> > l) 

= «5( W '),5(a;)«5( W - W ')» 

= (%),%)) = (l) s - (11) 
Now it follows that the expected value of the frequency of a signal S (u) can be written as 

M = ,«,(n>) a(f) . 

From this result, it follows that 



R>= .(^(vr) , (12) 

which can be proved by induction. If g (t) is an analytical function, it follows that 

( 9 M>= , ( ,,( 9 (w)) i(t) . (13) 

Thus, to calculate the average frequency of a function, we do not have to calculate the Fourier 
transform. Rather one simply calculates derivatives of a function and then integrates. 
The frequency translation operator has exactly the same effect: 

e i9f S (w) = S(u + 9). (14) 

For a complex signal, s(t) — A (t) e"^, 



SO 



e* w *(t) = ( 0' (t) - i-^ ) s(t) (15) 
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since the second term in the integral is a perfect differential. The average frequency is the 
derivative of the phase, over the density over all time. Thus the phase at each time 

must be instantaneous in some sense, i.e. Ui (t), so we can make the identification that 
uii {t) = (t). Similarly, we can show that 



S{w) 



(17) 



The covariance of a signal might be thought of as the "average time" multiplied by the 
instantaneous frequency or (t$ (t)) s = (t$' (t)) A ^ .When time and frequency are uncorre- 
cted with each other, then it is reasonable to expect that {M' (t)) = (t) (u>), so the difference 
between the two is a measure of how time is correlated to the instantaneous frequency. Thus, 
the covariance of the signal is 



Cov tu = (t# (t))-(t) (u), 



(18) 



while the correlation coefficient, r, is r = Covt ^ which is the normalized covariance. Real 
signals have zero correlation coefficients as do signals of the form A (t) e twot or S (u) = 
A(u>) e tuJt °, so signals with complicated phase modulation have a non-zero correlation coef- 
ficient. 

When dealing with more than one operator acting on a signal, we must be able to interpret 
the action of multiple operators such as AB acting upon signals. Here AB is taken to mean 
A acts on the signal followed by B acting on the signal. The commutator of A and B is 



A,B 



= AB - BA. 



(19) 



For example, the action of the time and frequency commutator on a signal is 



f , VV] s (t) = (fVV - VVT ) s (t) = is (t) . 



(20) 



This is analogous to the same result in quantum mechanics where the commutator of the 
position and momentum operator is equal to i when h — 1. The scale operator C is defined 

as 



<"5 



1 / d d 

2 V It It 



(21) 



It can also be written as 



c = rw + 



(22) 



C has the property that it transforms a signal s (t) according to 



(23) 



for a scaling parameter a. Thus, the wideband ambiguity function can be written as 

XwM^/e^) , (24) 

\ / s(t) 

the expected value of the scaling and translation operators for a signal s(t)e~ ltn f , which is 
equivalent to maximizing the signal to noise ratio (SNR) at the receiver. We explore what 
physical interactions, expressed in terms of operators, can be maximized. 



III. PHYSICAL INTERACTIONS 



While the primary scatterer produces the usual Doppler velocity and delay which is equiv- 
alent to the range, the operator viewpoint may hold some promise for finding interactions 
between the radar signal and the target that extend beyond considerations of position and 
velocity related criteria. Additional scatters can induce secondary characteristics into the 
return signal, such as micro-Doppler, which can be incorporated into the design of a receiver 
to maximize the possibility for detecting these types of secondary target induced charac- 
teristics. In addition to a scalar signal, higher dimensional waveform interactions can be 
considered as well, such as how the polarization of materials affects the waveform. The cross 
ambiguity function (CFA) symmetric form is defined as 



oc 



Xr >,r)= / e -^q*(t+-)s(t--)dt, (25) 



— oo 



where s(t) is the transmitted signal, while q(t) is the correlation signal and r is the delay 
parameter. This is the traditional form for the CFA. Instead of this form, a new type of 
CFA is proposed based on quantum mechanics. 

Any signal can be expressed as a complex vector. A new approach to signal amplification 
is presented here based on work by Aharonov on amplification of the measurement of some 
operators in quantum phenomena l]. Since any quantity that involves the usage of expected 
values of complex signals can be expressed in the same mathematical form as the quantum 
mechanical approach to signal amplification, the Aharonov approach suggests a potential 
candidate for the signal amplification that is similar to a CFA. The classical equivalent to 



this is what we choose to call cross correlation signal amplification. The definition of the 
cross correlation amplification of an observable A by the waveforms and |\I>/) is: 

. v (*f\A\*i) 

Across )i ~ (26) 

where both |\J/j) and f) are normalized. Now, the obvious question is how does the cross 
correlation measurement of an observable / (^A cross ^ differ from that of a normal observable 
i? 

Note that < (*/|*/) = 1, by the Cauchy-Schwartz inequality, so 

\(* f m\ <1. Thus, 

1 

and the effect of the denominator is to " magnify" the numerator provided there is no counter 
balancing effect. Note that if A \ = X^l^i), so 

f [Across).- -A A , 

so there is no effect. When there is not this cancellation effect, there can be a magnification, 
in some sense of the measurement of an operator. For an electromagnetic wave, the operator 
interactions can be treated as either two by two or four by four matrices. We consider only 
the two dimensional case. 



A. Multi-dimensional Interaction Operators 



Thus, the signal can be assumed to be of the form: 

l*i(*)> = 

and the cross correlation signal is: 

I*/ (*)> = 



E\ (t) 

EHt) 



(27) 



E{ (t) 
E{ it) 



(28) 



where the £"s can be real or complex. An interaction with a scattering object can be thought 
as a matrix, Ms, which acts on (t)) to give a return signal (t)), so 



\V R (t)) = M s \% (*)) 



(29) 



(weak) 



The cross correlation measurement amplification of operator M s is 

iKA v _ {9f(t)\* R (t)) 

f{Mcrossh ^ f{t)l%{t) y 

This example of amplification, which is analogous to spin systems in quantum mechanics, 
applies to polarimetric radars. Consider the four polarization matrices: 



00 



a 2 = 



1 




1 




, Ox = 




1 




-1 



1 




-i 








1 




i 



(30) 



(31) 



The first operator, ctq, acting on |\Pj(t)) is the identity, so it is equivalent to the previous 

2 ■ 2 

no amplification case. Now, if the waveforms are normalized, \E\ (t)\ + \E\ (t)\ = 1 and 
E[(t) \ E[(t) 2 = 1, so 



tan 9 



E\{t) 



and tan9' 



E{ (t) 



thus, we have 



E\ (t) E{ (t) ' 

nl E\ (t) E{ (t) 
tan^tan^'= - 2 



E\ (t) E{ it) ' 

Note, that we treated amplitudes as real so the angles are real, this is not necessary since 
complex angles are possible. The introduction of a complex angle as well would introduce 
a second term which is imaginary that would produce an additional effect on the imaginary 
component only. This possibility will be discussed in a future paper. 
Now, 

E\(t) 



°i l*< (*)> 



so 



(g, (t)!^ (*)) 

<*/ (*)> 



-E\ it) 

E{ (t) E\ (t) - E\ (t) Ej (t) 
E[ (t) E\ (t) + E\ (t) E{ (t) 
1- tan 9 tan 9' 



1 + tan 9 tan 9 



(32) 



When 9 ->■ -\, 



(VfWlhlVijt)) 1 + tan^ 
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so there can be amplification. In addition, we have 



* 2 |*i (t)) = 



E\{t) 
E\(t) 



so 



(* f (t)\o 2 \% (t)) _ E( (t) g (t) + E\ (t) E{ (t) 

(*/ (*) l*i (*)> ~ (t) ^ (t) + ^ (t) E{ (t) 
tan 0' + tan 9 
~ (tan ^ tan 0' + 1) 

sin (0 + 9') 
~ cos {9 - 9') ' 



When 0-0' ->■ f, 



<*,(*)! *■,!*,(*)> _ sin (26K + e + f ) 



e->0 C os( £ +f) 



— > oo, 



so amplification is possible. Finally, we have 



a 3 \^ l (t)) =i 



-E\ it) 
E\ {t) 



so 



(t)\ a 3 |* 4 (t)> -tff (t) (t) + ^ (t) (t) 



<*/(*) l*i(*)> 



£f (t) E\ (t) + ^ (t) E[ (t) 
itm (0-9'), 



tan a ± tan /? 



(33) 



(34) 



by using the trigonometric identity 

tan (a ± P) , 

1 =)= tan a tan p 

So amplification occurs as (9 — 9') — >■ |. Thus, the non-trivial operators <5"i, <t 2 , <7 3 can be 
amplified under the right conditions for the components of cross-selection waveforms. 
There are four additional operators to consider: 



li 



Pi 



12 



1 


1 
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21 



and 



Now, 



so 



22 





1 


1 



Ail*, (*)> 



<*/ (*)l*i(*)> 



Ej(t) 


E[ (t) E\ (t) 



E[ (t) E\ (t) + E\ (t) E{ it) 
1 



(35) 



tan 9 tan 9' + 1 ' 

Note, the denominator goes to zero as tan 9 tan — > — 1, while the numerator remains finite, 
so amplification is possible for this operator. Also, if Ai is replaced by a constant aP n , the 
amplification effect works as well. Now, 



Pl2\*i (*)> 



SO 



Pl2|*i(t)) 
<*/ (t)\%(t)) 



Ei(t) 



^ (*) E\ it) 



E( (t) E\ {t) + E\ it) E{ it) 
tan 9 



. - (36) 
tan 9 tan 0' + 1 

Note, the denominator goes to zero as tan tan 0' — >■ —1, while the numerator remains finite, 
so amplification is possible for this operator. Also, if P 12 is replaced by a constant aPi 2 , the 
amplification effect works as well. Now, 



(*)) = 



so 



<*/ (*)l*i(*)> 





£j(t) 



4 (t) £j (*) 



flf (t) (t) + (t) ^ (t) 

tan 0' 



tan tan 0' + 1 ' 



(37) 
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Note, the denominator goes to zero as tan 9 tan 9' — > — 1, while the numerator remains finite, 
so amplification is possible for this operator. Also, if P21 is replaced by a constant (2P21, the 
amplification effect works as well. Now, 

P22\% (*)) = 

SO 

P 2 2|*t(t)) _ E{(t)EW) 
(*/ (*) l*i (0) ' Pf (t) E\ (t) + E\ (t) E{ (t) 

tan 9 tan 9' 

~ tan 6> tan 9' + 1 ' ^ ^ 

Note, the denominator goes to zero as tan 9 tan 9' — > —1, while the numerator remains 
finite, so amplification is possible for this operator. Also, if P22 is replaced by ap22, the 
amplification effect works as well. 

Note, we have provided the necessary conditions under which these operators can be 
amplified, but they are not sufficient. Sufficiency comes when waveforms can be shown to 
obey the conditions the angles obey to produce amplification. These conditions must be 
shown to be satisfied by specific waveforms or classes of waveforms. In addition, noise has 
to be brought into the mix. 







B. Scattering Operators 

The scattering operators for five specific structures are examined from the viewpoint 
of amplification of operators. These scattering operators special cases, two dimensional 
matrices, of the more g eoe r a, operate, fcd—^atrtces, found in Conettfl. 

1. For a sphere, a plane, or triangular corner reflector oriented horizontally, the scattering 
matrix is: 

= <t . (39) 

Since this is the identity, there is no amplification effect. For a sphere, a plane, or 
triangular corner reflector vertically polarized, the scattering matrix is: 

= i**, (40) 



S(h,r) 



1 
1 



S (v, r) 



% 

1 



12 



so it can be amplified. (Note h stands for horizontal polarization and v stands for 
vertical polarization.) 



2. For a dipole oriented along the vertical axis is : 



S(h,r) 



1 




11) 



and 



S(v,r) 



1 -i 
-i 1 



— <J2 + -(To. 



S (h, r) can be amplified, while the first term of S (v, r) can be amplified. 
3. For a dipole oriented at the angle a from the positive horizontal axis: 



S(h,r) 



cos 2 a | sin 2a 
| sin 2a sin 2 a 



- sin lab 2 + cos 2 aPu + sin 2 aP 2 2, 



and 



S (v,r) = - 



e i2a _• 



—i e 



-ila. 



—a 2 + e l2a Pn + e- l2a P 22 . 



(41) 



(42) 



(43) 



(44) 



Since a 2 and aPu + j3P 22 ^ cr can be individually amplified, then 



aP n + PP- 



22 — — 



a 

p 



where S = S^which implies that aPu + pP 22 is Hermitian. 



4. For a dihedral corner reflector oriented along the horizontal axis: 



S(h,r) = 

and 

S(v,r) = 



1 

-1 

1 
1 



a . 



The first operator can be amplified and the second can't. 



(45) 



(46) 



(47) 
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5. For a right helix oriented at an angle a from the positive horizontal axis: 



S(h,r) 



-i2a 



1 -i 
-i 1 



-i2a 



[a - ia 2 \ 



(48) 



and 



S (v, r) 





e~ l2a 



e -i2ap 



22- 



(49) 



Clearly the matrix S (v, r) is amplified. For S (h, r), although a® is not amplified, the 
component e ^ [^2] is amplified relative to it. 

6. For a left helix oriented at an angle a from the positive horizontal axis: 



-Via 



S{h,v) 



1 i 
% -1 



-i2a 



-i2a 



e „ ze 

-TT~ a l H S — ^2- 



(50) 



Clearly this operator can be amplified. 



C. Single Dimensional Interactions with Signals 

The goal of receiver design is to maximize the response of a receiver with respect to the 
return signal sr (t). The functional form is sr (t) = s(at + r) where r is the (delay) time it 
takes the signal to reach the target and return to the receiver, and a is dilation of the time 
axis due to the motion of the object. This is accomplished by taking the inner product of 
sr (t) with s* (t) and integrating, so we are computing the Fourier transform of the product 
s*(t)s(at + r): 



A 



A' 



( s(t)e Uu) , s(at ±t)) = ^=( e iaC e trW 

V7T \/TT \ 



(51) 

which is the expected value of the operators for scale e and the operator for time shift 
girw Tryi n g t maximize the reception SNR has led to the ambiguity function which can 
be interpreted as the expected value of two specific operators for a given signal s(t). 

The non- uniform Doppler effect can be used to illustrate this operator viewpoint. The 
effect of non-uniform Doppler on the radar waveform can be determined by the application 
of the relativistic boundary conditions to the D'Alembert solution to the wave equation [ 
The scattered waveform in terms of the incident waveform becomes 

2r(r) 



g(j) - f [r 



(52) 
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For a dynamic system characterized by single parameter a, then a dynamic variable u evolves 
along a path in configuration space. The configuration of the system describes a curve along 
a. Consider the commutator equation 



du 
da 



u, G 



(53) 



Here, G generates the trajectory u = u(a) and a can be viewed as geometrical parameter. 
Expanding u(a) in a Taylor series yields jl7j| a Taylor series, thus the generator equation can 
be used to replace the dynamics with the operator equation 

,,2 rr 

exp (otG) u(a), a „ . 



U\Oi) 



Uq + a 



u,G 



a 

a=o + 2! L 



u, G 



G 



+ 



a=0 



(54) 



For physical systems, it is evident that the generator of dynamics is time, so any function 
of time can be thought of as being generated by an operator, G, acting on u (t), so it can be 
thought of being "generated" by that operator. It is evident how to "generate" any function 
of a parameter using operator methods 7]. For a given r(r), we can assume it is generated 
by a equation such as r(r) = exp [kGj x(r)\ T=0 , so 



f r 



2r(r) 



fir - a exp G r(r) | T=0 ) = exp [oH )s (r) , 



(55) 



where H depends on the specifics of the interaction. For example, H would be a comb 
operator in the frequency domain for a periodic function. In this case, we are estimating the 
expected value ^exp ( r ^) ^ a ^ the receiver. Since any scalar interaction on the waveform 
can be thought of as the action of an operator on the broadcast waveform, a more general 
ambiguity function can always be defined as 

Xo(^,t) = ^L f e~ ituJ s*(t)exp(rn) s(t)dt = (e~ m s*(t),exp (rn) s(t)) . (56) 



For the remainder of the discussion, we assume the signal is not normalized. The typical 
signal processing application is to minimize the effect of the noise h so as to maximize the 
signal-to-noise ratio (SNR) for a received signal y. In order to understand how to do this, 
one uses a linear model for the combination of signal plus noise y — s (t) + n. The response 
to an input f(t) of a system function h(t) is the response g{t), which at a time to is 



g(tc 



s(t)h(t - t )dt 



1 

2^ 



S{u)e wt0 H(u)du. 
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Here we wish to determine the maximum value of g(to)] this allows us to maximize SNR 
depending on which of several integral constraints that are specific to the problem being 
considered. The SNR depends on the mean squared constraint under consideration: it could 
be based on the energy spectrum ^(u;)! 2 , it could be based on the constrained energy 
spectrum \S (u)\ 2 \R(u)\ 2 , it could be based on multiple constraints such as higher order 
moments of the energy spectrum, or it could be based on amplitude constraints. Each 
constraint leads to a different choice for the system response function H (u>). 
If we have a specified energy 

1 f°° 

E = (s(t),s(t)) = —J IS^fdu, (57) 
then by the Cauchy-Schwartz inequality 

/oo 2 /»oo poo 

S(u)e iojto H(uj)duj < \S(u)\ 2 du \e wta H (u)\ 2 du 
-oo J —oo J —oo 

oo 



<— I \e tuJt0 H(u)\ 2 du. (58) 

27T J. 



oo 



The inequality becomes an equality only if 

S{u) = ke iut0 H*(u), (59) 

so the maximum value for g(t ) is obtained by the choice 

s(t) = kh*(t -t) (60) 

since H* (oo) i — > h*(—t) and k is an arbitrary constant. For a linear system y(t) — s (t) + 
n(t) with an impulse response h(t), the output u(t) is 

u(t) =y(t)*h(t) =u s (t)+u n (t) (61) 



where u s (t) — s(t) * h(t) and u n (t) — u(t) * h(t). Now the response to the signal s (t) is 

f oo 

2tt 



u s (t ) = / S(u) e^H (w) du, (62) 



so 

\u. 



/OO poo 
S (u) e lujt " H (u) dco / S (u) e iult » H (u) du; (63) 
-oo J — oo 
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this is what Papoulis has called the Matched Filter Principle jl^. From the operator per- 
spective, the operator acting on the signal should replace the operator acting on the system 
response function in this argument, so 

g(t )= [ s(r)exp(aH) h(r -T )dr = — [ S (u) e i0JT ° R (u) du) . 
J-oc v J 2tt J_ oa 

where 

(V^exp (all) h(r - r )^ = (e lbJT ,exp (ail) e r °^/i(r)^ = e iulT0 (e iUT ,exp (ait) h(r)) = e iulT °R{uj) 

(64) 

since the operators commute. When 

exp (aii) = e iad e iT ™, (65) 



and the optimum choice is a rescaled version of the transmitted signal time scale to — > at±r, 
the wideband matched filter. 

The Matched Filter Principle is quite general and can be used to introduce a variety of 
constraints, which are equivalent to a cost function minimization approach. For example, if 
one wanted to maximize the response to the derivative of the energy E\, while requiring the 
energy to be normalized, then one has 

/OO -| poo 

\ s ' (r)| 2 dr = — u 2 \F{u)\ 2 du (66) 
-CO ^ 7I " J — OO 

so \s' (t)| < \[E[ with equality at time to if s (r) = \[E[ exp O^rj'j- In general, using this 
approach, arbitrary constraints can be considered. If we have a signal exp (ail) s (r) of 
where the energy of s (r) is E, that we want maximize the system response g (to) of the 
system h (t), then to obtain the maximum subject to the constraints 

s(t)$i(t)dt = #i, (67) 

where the functions (t) and constraints i?j are given. Then, with the definition 

Ui{t) = ^ {t) ~wr (68) 

that the constraint equation becomes: 

o 

s (t) Ui (t) dt = (69) 
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because S(0) is the area of s (t). Thus, it follows that the system response is 



/CO . . n 

s(r) exp (aVA e T °^h{r) + ^/3 i M i (r 
-°° V i=i 



for arbitrary Therefore, <?(ro) can be bounded by 



\9(roW< 



Equality is achieved if 



/oo n 
exp (aVA e T °^h{r) + 
-°° L i=l 



S T 



n 

exp («7i) e^/^+J^* (r) 

i=l 



dr, 



dr. 



(70) 



(71) 



(72) 



This gives a method for choosing the correlation waveform to achieve maximum response 
for a given set of constraints. 

IV. CONCLUSIONS 



The operator method is a much richer way to look at the radar measurement problem 
because of its ability to produce a wide variety of distributions associated with the infor- 
mation contained in a signal. In particular, it is possible to put the ambiguity function in a 
wider context as part of a general theory of measurement. There is a much greater freedom 
of description of the same physical situation which suggests that we can find information 
present in waveforms that a waveform designer would not think to look for. This approach 
to incorporating quantum mechanical ideas has been championed by Baraniukj^] js] recently 
by extending the Hermitian operator approach in quantum mechanics to unitary operators 
in signal processing. The specifics of the type of operators matter relative to the physics 
of the interaction of the target with the waveform, so this may be important for future 
extensions of this work. 
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